Enhancement of thermal entanglement in two-qubit XY models 
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We analyze conditions leading to enhancement of thermal entanglement in two-qubit XY models. 
The effect of including cross-product terms, besides the standard XY exchange interactions, in the 
presence of an external magnetic field, is investigated. We show that entanglement can be yield at 
elevated temperatures by tuning the orientation of the external magnetic field. The details of the 
intrinsic exchange interactions determine the optimal orientation. 
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I. INTRODUCTION 

Entanglement implies strong correlations between spa- 
tially separated quantum systems that can not occur clas- 
sically. By virtue of its nonclassical nature, entanglement 
has started to be seen in recent years as a phenomenon 
at the core of the future technology, namely, as a poten- 
tially useful resource for quantum information process- 
ing. Nevertheless, beyond this exciting possibility, quan- 
tum entanglement deserves further investigation as it is 
not a rare phenomenon but it is generic for the states of 
interacting many particle systems and, therefore, ubiqui- 
tous. 

In the preparation of entangled states, the external 
control performed through changes in temperature, elec- 
tromagnetic fields, etc., should not introduce significant 
levels of decoherence and, at the same time, the required 
values of the external control parameters must fall within 
realistic working ranges. Due to the usually large de- 
coherence times, the spin is a property that has been 
exploited in many implementations requiring entangled 
states [1]. As a consequence, most of the solid state 
proposals lie upon two-particle exchange interactions of 
the Heisenberg type [2], so that, the intrinsic two-body 
Hamiltonian operator is of the form 



Hi. 



E 



(1) 



,(2) 



(1) 



where a^"^ are the Pauli matrices of qubit n and J* the 
coupling coefficients {i = x,y,z). For instance, one finds 
the isotropic Heisenberg model (J^ = = J^) in spin- 
coupled quantum dots [3] and donor-atom nuclear spins 
[4]. Its XY version [J^.J^ ^ = 0) is found in quan- 
tum dots in cavities [5] , excitons in coupled quantum dots 
[6], atoms in cavities [7] and quantum-Hall systems [8]. 



Ising-like systems (J^ 7^ = = 0) have also been 
considered [9] . 

The entanglement of two-qubit systems for the XY 
model, in the presence of a constant external magnetic 
field BzZ, has been investigated recently [10,11]. It was 
shown that the thermal entanglement of anisotropic sam- 
ples can be manipulated through changes in Bz ■ It is our 
purpose here to analyze the effect that the inclusion of 
cross-product terms of the form crp-' (g) (Tj^\ with i 7^ j, 
yields upon thermal entanglement. These terms model 
weak ferromagnetism and originate from the spin-orbit 
coupling [2,12]. Also in the scenario of fermion quan- 
tum circuits, the anisotropy in the exchange interaction 
may be accompanied by cross-product terms [13]. For 
external control, we will introduce a constant external 
magnetic field B = {Bx,By, B^)- 

The degree of entanglement will be gauged through the 
concurrence [14], a measure that, for two-qubit systems, 
gives the exact frontier between separable and entangled 
domains. Briefiy, the concurrence C is calculated as 



C 



:{2'max{\i\ — '^^Xi, 0} 



(2) 



with {Xi} the square roots of the eigenvalues of the 
matrix R = pSp*S, where p is the density matrix. 



S = o'y"^'' ® ay"' and stands for complex conjugate. A 
nonzero concurrence means that the two qubits are en- 
tangled, in particular unitary concurrence corresponds to 
maximally entangled states. 

The paper is organized as follows. In Section II we 
describe the system. The results are presented and dis- 
cussed in Section III. Finally, Section IV contains the 
concluding remarks. 
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II. THE SYSTEM 

In order to study the effect of extra cross-product terms 
in the XY model, let us consider the following internal 
Hamiltonian 

Ho=\Y. <^'o^®-?, (3) 

with a^^ = (1 + l)J, = (1 - 7) J, where J G 3? is 
the coupling constant and 7 e [—1, 1] is an anisotropic 
parameter. Also, let us define the cross couplings as 
a^y = (1 + 7')-?!: and a?'^ = (1 - i)K with if G 3? 
the cross coupling constant and 7' G [—1,1]. Then, Ho, 
written in terms of the raising and lowering operators 
= (ai")±i4"))/2, reads 

Ho = (7 J - iK) u^^^ ® + ( J + iiK) a^^^ ^ a^^^ 

+ (J - i^K) ® (7^'^ + (7 J + iK) ® . (4) 

It is clear now that the inclusion of cross-product terms 
is equivalent to setting complex coupling constants. 

In the presence of the magnetic field B, the complete 
dimensionless Hamiltonian becomes 

H = Ho + l ^ ^.(af'eaf)). (5) 

i=x,y,z 

We consider states of equilibrium at temperature T. 
Thus, the density matrix is given by p = exp{—H/T)/Z 
with Z = Tr{exp{—H/T)) where we have already set the 
Boltzmann constant fc^ = 1. The entanglement at these 
states is called thermal entanglement [15]. 

III. RESULTS 

Let us analyze first the case B± = 0, B\^ = 0, 

where is the component of the field in the a;y-plane. 
The eigenstates of the density matrix p are 

= [{J + ^l'K) lU) ±^\l^)]/N^ (6) 

(with eigenvalues ±/x, respectively, being jj = 
V + (7'^)^ ) and 

I <t>t) = [ilJ-iK) ITT) + {±X-B^) \U)]/Nx (7) 

(with eigenvalues ±A, being A = yV^ + fiJ, = 
{jJ)'^+K'^ ), where A^^^, N\ are normalization constants. 
The eigenvalues of matrix R are given by 

A?,2 = exp(±2M/T)/z2 

A^_4 = (1 + 2a;2 ± 2xVl + ^)/Z^ (8) 

with X = {v/X) sinh(A/r) and Z = 2cosh(i//r) + 
2cosh(A/T). Then, the concurrence C is straightfor- 
wardly computed through Eq. (2). At T ^ 00 or 



B\_ 00, the entanglement vanishes (C — > 0). At T = 0, 
the concurrence is 

f 1 M > A, 

C(r = 0) = <^ (l-i//A)/2 for /x = A, (9) 
[ u/X /i < A. 

Hence, for T = 0, there is a jump in the concurrence at 
/U = A similarly as already observed in the literature for 
the case K = 0,7 7^ [11]. For /x > A, the ground 
state is | which is fully entangled. For /i < A, 

the ground state, | 0~)) is partially entangled unless 
B± = in which case it becomes fully entangled. But, 
at ^ = A, the ground state disentangles partially(totally) 
for BjL + 0(=0). 

The results are not dependent on the sign of the cou- 
pling parameters, therefore, in particular, they are valid 
both for the antiferromagnetic and ferromagnetic cases. 
Also they are independent of the sign of B^. 

The picture in the absence of magnetic field {\B\ = 0) 
is recovered by setting B± = in expressions (6)- (9). 
In Fig. 1 we show T*, the temperature at which the 
concurrence vanishes, as a function of parameters 7', for 
different sets of {J,K), in the absence of any external 
magnetic field. The concurrence is null above T*. If 
7'^ > 1 - {J/KY, then, the threshold temperature T* 
grows with increasing 7'. For \B\ ~ 0, by interchang- 
ing the values of parameters (if, 7') with those of (J, 7), 
one obtains the same results. Taking this into account, 
notice that, even for the Ising case (I7I = 1), there can 
be entanglement at finite temperature in the absence of 
magnetic field {\B\ = 0), as soon as K is large enough. 

Fig. 2 exhibits the effect of the external field B± on an 
isotropic sample (7 = 0). By increasing |i?_L|, it is pos- 
sible to increase the threshold temperature T* provided 
K ^ 0. In reference [11] it was shown that anisotropic 
exchange interactions are required for allowing control of 
T* by varying the intensity of B± . However, the presence 
of cross-terms enables control of the threshold tempera- 
ture by varying the intensity of the magnetic field in the 
^-direction, even in the absence of anisotropy. 

If i3|[ = 0, the subspaces spanned by {|tT)i lii)} and 
by {|Ti))|Ti)} remain uncoupled, while for i^y 7^ 0, they 
mix. Numerical calculations are performed in this case. 
As expected, the results are only dependent on the to- 
tal component of the magnetic field in the xy-pl&ne, B|| . 
Plots of the concurrence as a function of T for different 
values of By , being B± = 0, are presented in Fig. 3. a. For 
T = 0, there is a jump in the concurrence at B|| = \/2J 
when 7, if = 0. In this case, one obtains 

r 1 By < V2J, 

C{T = 0) = I 2/3 for B,[ = V2J, (10) 
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with D = (^1 + (2B||/J)^ - 1)/(2B||/J). 

For B = B±z and = 0, it was already shown in 
the hterature [10] that, when B± > J, although there is 
no entanglement (C = 0) at T = (also in agreement 
with Eq. (9)), the concurrence presents a maximum at 
finite T. In fact, at T = 0, the ground state is the dis- 
entangled state I J, J,), that mixes with the maximally en- 
tangled ones as T increases. A similar effect is observed 
for B = S|| and 7, = (see Fig. 3. a), although in 
this case C(T = 0) 7^ consistently with Eq. (10). For 
By < V2J, the concurrence C decreases monotonously 
with T, whereas, for S|| > \/2J, it presents a local max- 
imum at finite T. 

A plot of C vs. \B\ is shown in Fig. 3.b for different 
orientations, at T ~ and 7, ii' = 0. For all orientations 
there is a jump at T = 0. It occurs when 

2Bi + B| = 2j2. (11) 

Differently from the case B = B±z, if B\\ ^ 0, then, the 
concurrence falls to a non null value above the singular- 
ity. Thereafter, C tends smoothly to zero with increasing 
intensity of the magnetic field. 

In Fig. 4 we exhibit T* as a function of the intensity 
\B\ for different orientations of the field when 7 = and 
K = {). Notice that if the magnetic field is oriented in 
the ^-direction, then the threshold temperature does not 
depend on the intensity of the field. However, when the 
field has a non null component in the xy-plane (By ^ 0), 
then, T* increases for sufficiently large B\^. By compar- 
ing Figs. (2) and (4) notice that B\^ ^ yields an effect 
qualitatively similar to that due to Sj_ when the system 
has cross-terms (K 7^ 0). A similar effect is also ob- 
served with B^ when the system is anisotropic (7 ^ 0) 
[11]. But, even in the absence of anisotropy or cross- 
terms the external magnetic field can be used to control 
the entangled domain by switching on the || -component. 

IV. CONCLUSIONS 

To summarize, we have investigated the thermal en- 
tanglement of two-qubits XY systems with additional ex- 
change cross-terms (K ^ 0) as well as with anisotropic 
couplings (7 7^ 0). Arbitrary orientations of an external 
magnetic field B were considered. We have found that, 
even in the absence of anisotropy, cross-terms enables 
control of the threshold temperature T*, which defines 
the domain of entangled states, by varying the intensity 
of the magnetic field in the 2-direction. Moreover, even in 
the absence of anisotropy or cross couplings, the degree 
of entanglement can be manipulated by suitably chang- 
ing the relative orientation of the sample with respect 
to a constant magnetic field. As a perspective, it would 
be interesting to investigate the effect of increasing the 
number of spins. 



Since thermal entanglement is a natural type of en- 
tanglement for a system embedded in a thermal environ- 
ment, we expect that the present results could be useful 
for solid state applications. In particular wc expect that 
the understanding of thermal entanglement could be ex- 
ploited within the context of quantum communication 
and information processing, e.g., along the lines of quan- 
tum heat machines [16]. 
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FIGURE CAPTIONS 



Figure 1: Temperature T* above which the concurrence 

vanishes as a function of 7', in the absence of magnetic 
field = 0, for 7 = and different values of {J,K) 
indicated in the figure. 

Figure 2: Temperature T* as a function of B±, for the 

parameters indicated in the figure. The entangled do- 
main is the low-temperature region bounded by the curve 
and abscises axis. 

Figure 3: Concurrence C as a function of T (a) for dif- 
ferent values of _B|| and B± — and C as a function of 
|_B| at T ~ (b) for different orientations of the magnetic 
field. B± = \B\ cos 6* and = sinO, where the values 
of angle 9 are indicated in the figure. Here 7 = 0, ii' = 0. 

Figure 4: Temperature T* vs \B\ for the orientations of 

the field indicated in the figure. The entangled domain 
is the low-temperature one. Here j = 0, K = 0. 
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